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/—* In this work the junction conditions between the exterior Reissner-Nordstrom- Vaidya space-time 

with the interior quasi-spherical Szekeres space-time have been studied for analyzing gravitational 

^vj ■ collapse in the presence of a magneto-hydrodynamic fluid undergoing dissipation in the form of heat 

flow. We have discussed about the apparent horizon and have evaluated the time difference between 

^ , the formation of apparent horizon and central singularity. 
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I. INTRODUCTION 

O 

In Einstein gravity, gravitational collapse is an important problem for astrophysical objects. Usually, the 
K' formation of compact stellar objects such as white dawrf and neutron star are preceded by a period of collapse. 
Hence, for astrophysical collapse, it is necessary to describe the appropriate geometry of interior and exterior 



regions and to determine proper junction conditions which allow the matching of these regions. 



The study of gravitational collapse was started long ago in 1939 by Oppenheimer and Snyder [1] to solve an 
idealized problem of gravitational collapse preceded by assuming a spherical symmetric distribution, adiabatic 
flow and the equation of state in the form of dust. Misner and Sharp [2] have introduced pressure gradient 

^-H ■ terms in the equation of motion. Outgoing radiation of the collapsing body has been considered by Vaidya 
[3], Santos [4] and collaborators [5 - 10] and the correspondence to adiabatic distribution of matter has been 
considered by Misner [11] and collaborators [12, 13]. Models of radiating collapse with viscous fluid were 

r^. ' presented by Lake [14] and Santos [4]. Ghosh and Deshkar [15] have considered collapse of a radiating star with 
a plane symmetric boundary (which has a close resemblance with spherical symmetry [16]) and have concluded 
L| with some general remarks. 
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Most of the studies has considered a collapsing star whose interior geometry is spherical. But in the real 
astrophysical situation, the geometry of the interior of a star may be quasi-spherical in form. Previously we 
have discussed [17] junction condition between interior Szekeres model and exterior plane symmetric Vaidya 
space-time. Oliveira et al [9] have proposed to study the junction conditions of a collapsing non-adiabatic 
charged body producing radiation. The interior space-time V~ was modeled by a magneto-hydrodynamic 
fluid undergoing dissipation in the form of a radial heat flow and the exterior space-time V + was described by 
the Reissner-Nordstrom- Vaidya metric which presents a spherically symmetric electric field and a radial flow 
of unpolarized radiation. Now in this paper we will consider the interior space-time V~ is quasi-spherical 
Szekeres model and the exterior space-time V + is Reissner-Nordstrom- Vaidya metric in the presence of 
electro-magnetic field. 
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II. JUNCTION CONDITIONS 

Let S be time-like 3D hypersurface which separates two AD manifolds V~ and V + . The interior 
AD manifold V~ is given by the quasi-spherical Szekeres space-time with line element 

V~ : ds 2 _ = -dt 2 + e 2a dr 2 + e 2f3 {dx 2 + dy 2 ) (1) 

where a and (3 are functions of all space-time co-ordinate variables. 
Now for the transformation 



X = cot — cos0, y = cot — sin0, 



the metric (1) in V becomes 



dsl = -dt 2 + e 2a dr 2 + -e 2fl cosec-{d0 2 + sin 2 9 dej) 2 ) (2) 

We assume that the source for the Einstein's field equations in V~ is given by 

Gfj, v = T^„ + E^ v (3) 

The stress-energy tensor of a non- viscous heat conducting fluid has the expression 

T^ u = (p+p) u^Uv+p g^ v +q^u v + q v u^ (4) 

where p, p, (/„ are the fluid density, isotropic pressure and heat flow vector respectively and u^ is 
the four velocity. We take the heat flow vector q^ to be orthogonal to the velocity vector u^ i.e., 
q^ u^ = 0. For co-moving co-ordinate system, we choose u^ = (1,0,0,0) and g M = (0,q r ,q x ,q y ) with 
q t =q t (t,r,x,y), i=(r,x,y). 

The energy- momentum tensor E^ v of the electro-magnetic field is given by 

1 



Env = — 

p Air 



F v F a F vS F x 



(5) 
where i* 1 ™ is the electro-magnetic field tensor satisfying Maxwell's equations 



^ = 0i/, M - </y,i/ (6) 

and 

Ff„ = 4tt J^ (7) 

where cj)^ is the four-potential and J M is the four-current. 

Since the charge is at rest in this system, so there will be no magnetic field in this system. Thus, we may 
choose the four-potential and four-current as 

0„ = (0(t,r,s, W ), 0,0,0) (8) 

and 

J" = a u" (9) 



where a is the charge density. From equation (6) using (8) we have the only non-zero component of the 
field tensor as 



Foi=-F 10 = -^ (10) 



Also from equations (7) and (9) we obtain 

d 2 s 



gr > -«'i-W« (ii) 



and 

8 2 c 



d^=0 (12) 



drdt dr 

where ' ' ' and ' . ' denote partial derivatives with respect to r and t respectively. Integrating (12), we 
have 



and then from equation (11) one gets 



^- = K(r,x iy )e a (13) 



K' = 4ir<re a (14) 

For the interior space-time model (1) with matter field (4) and electromagnetic field (5), the Einstein's 
equation (3) can be explicitly written as the following : 

2 a $ + f3 2 - e- 2f3 (a 2 x + al + a xx + (3 XX + a yy + f3 yy ) + e - 2a (2 a' - 3 (3 2 - 2/3") = p - ±- K 2 (15) 

3 /? 2 + 2 (3 - e~ 2a I3 2 - (f3 xx + & m )e~^ = -p- ^-K 2 (16) 

a 2 + & + aj3 + (3 2 + + e- 2a (a'(3' - (i 12 - [3") - e- 2f} [a 2 x + a xx - a x f3 x ] - e~^a y p y = -p + —K 2 (17) 

a 2 + & + a$ + $ 2 + J3 + e- 2a (a'(3' - (3 12 - (3") - er 2 f 3 [a 2 + a yy - a y (3 y ] - e- 2l3 a x (3 x = -p + ^-K 2 (18) 

a x (-a y + P v ) + (3 x a y ~ a xy = (19) 

ot y (-a x + (3 X ) + [3 y a x - a xy = (20) 



a/3' + [5(3' + = X - q r e 2a (21) 



(d - /3)a x + u x +(3 x = q x e 2/3 (22) 

(a - $)a y + a y + $ y = q y e 2/3 (23) 



a x - X = 



(24) 



a y 0' - (3' y = 
where subscript stands for partial derivatives with respect to corresponding variables. 

For quasi-spherical case i.e., for (3 1 ^ , we have the solutions [17] 



(25) 



and 



e 13 = R(t,r) e v{r ' x ' v) 



R' + Rv' 



D(t,r) 
where R(t,r) and D(t,r) will satisfy the differential equations 



(26) 
(27) 



2 R R + (R 2 - D 2 ) + ( p + |- ) R 2 = f(r) 



and 



R D = f(r) e- 2a 
with f(r) an arbitrary function of r . The function v can be written in the form 



(28) 



(29) 



e~ v = A{r) (x 2 + y 2 ) + B^r) x + B 2 (r) y + C(r) 
where A,B\,B-2 and C are arbitrary functions of r along with the restriction 



(30) 



(B 2 + B 2 )-4AC = f(r) - 1 
From equations (21), (22) and (23) we have the components of heat flux vector as [17] 



(31) 



2 • 
Qr= R D ^ 



(32) 
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D 
D 



a r e 



(33) 



D -fi 

The exterior manifold V + is described by the Reissner-Nordstrom-Vaidya space-time having metric 



(34) 



V+ : ds 2 



2M(v) Q 2 



z z* 



2{jn2 i „• 2 n JJ.2\ 



dv z - 2dv dz + z 2 (d8 z + sin 



(35) 



This space-time represents a spherically symmetric field with an outgoing radial flux of unpolarized radiation 
and v is the retarded time. 



The surface of separation £ between these two manifolds V is characterized by r — r^ and the intrinsic 
metric on it (in the comoving co-ordinates of the interior space-time V~ ) is given by 

dsl = -dr 2 + R 2 ( T )(d6 2 + sin 2 d(f> 2 ) 
Now in the description of Santos, the Israel's junction conditions are 
(i) The continuity of the line-element 

(dsl)* = (ds 2 + )j; = ds 2 , (36) 

where ( )s means the value of ( ) on £ . 

(ii) The continuity of the extrinsic curvature over £ : 

[Kij] = [K±] - [Kr.] = 0, (37) 

where the expressions of extrinsic curvature in terms of the unit space-like normal vector to £ , n^ are 
given by [18], 






Kfj = -<^ - ntT%,-£-£- (38) 



Here £* refer to the intrinsic co-ordinates (r, x, y) on £ and x± > A* = 0, 1, 2, 3 stand for the co-ordinates 

in V r± . 

The boundary £ in terms of the co-ordinates of the interior space-time V~ is given by 

•F(r, t) = r — rs = (rs , a constant) (39) 

So one can write the normal vector n~ as 

n; = (0,e«,0,0). 
Now comparing the metric ansatz for £ and V~ for dr = , the continuity of the line element gives 

— = 1, TZ(t) = e li on r = r E (40) 

dr 

Further, the components of the extrinsic curvature for the interior space-time are 



K,. = 0, K aa = cosec 9 K, 



and K { , = 0, i ^ j 



\0 2 e 2 P~ a cosec^l) 



(41) 



On the other hand, in terms of the co-ordinates of the Rcissncr-Nordstrom-Vaidya space-time the surface 
£ can be characterized by 

F{z,v) = z-zv{v) = (42) 

Hence, as before the unit normal vector to £ is given by 



and the components of the extrinsic curvature for the exterior space-time are 



K+ 



v V z z z A j 



K+ = cosec 2 6 K+ = 



K?j=0, %*j 



■ i -i 2Mlv) . Q , 

V ( 1 z^- + ^ ) Z + ZZ 



(44) 



(here ' . ' stands for the differentiation with respect to ' r ' ). 

Now, the continuity of the metric ansatz for S and V + over z = zs(v) is given by 



ze(v) = R(t) 



) dz 

' dv 



2M(v) 



(£) 



v 2 )T. 



(45) 



Also the continuity of the components of the extrinsic curvature across S give 



p = q r e 



/(0 - 1 
R 2 



R 2 



sin 4 (9/2) e - 



R 4 ~ 8tt 



and 



,,. . i? Q 2 L> 2 i? i? i? 2 



(46) 



(47) 



which can be termed as the total energy entrapped inside the surface S . We note that on the boundary, the 
pressure can not be vanish in general. Thus for quasi-spherical shear-free distribution of a collapsing charged 
fluid undergoing dissipation in the form of heat flow, the isotropic pressure on the surface of discontinuity can 
not be zero. In particular, if fluid stops dissipation, q r = 0, the pressure can not vanish at the boundary, 
but in this case radiation can not exist and the exterior space-time V + is the Reissner-Nordstrom-Vaidya 
space-time. But in absence of isotropic pressure there may still be radiation on the boundary and the exterior 
space-time V + will still be Reissner-Nordstrom-Vaidya space-time. 

The total luminosity for an observer rest at infinity is [7] 



L r 



dm\ 
dv / , 



(48) 



If we now consider an observer on the boundary S then the luminosity for that observer is 



U 



dm 
dv 



(49) 



Now the boundary redshift of the radiation emitted by a star can be written as 



Zy 






dv 1 

At R + D 



- 1 



(50) 



Hence the luminosity measured by an observer at rest at infinity is reduced by the redshift in comparison to 
the luminosity observed on the surface of collapsing body. Also when R + D = then the boundary redshift 
attains unlimited value (i.e., Z-% — > oo ). 



III. APPARENT HORIZON AND ITS TIME OF FORMATION 

The first integral of equation (28) is given by 



R 2 = /(? 



F(r 



M " ± i^-pB^W 



R 24tt 

As D is regular and p blows up at the singularity [19], so let us assume 



(51) 



D = X(r)R, p 



m 

R 2 



The evolution equation (51) satisfies to 



R 



b'2 



F(r) K 2 R 2 A 



■R 2 



R 24tt 3 

Further, the time of formation of apparent horizon t a h(r) is given by 

R 2 (t ah (r),r) = l + f(r) 
Hence for equation (53) we have a cubic equation in R as 

T " S ^-"W- r ) - ^ + f(r))R(t ah (r), r) + F(r) = 



(52) 



(53) 



(54) 



(55) 



Thus apparent horizon exists if this cubic equation in R has at least one positive root. We shall show below 
different conditions for existence of positive roots. 



I. For 4j — -^ > and F(r) < — - ■'' 2 t , there are two positive roots of equation (55) and hence there 



3 24tt 

are two apparent horizons, given by 



Rc(r) = 2. 



(1 + /) 3 

(hi _ K 2 \ 
y 3 24tt/ 



-COS 



2 (r) \ 3 24tt 



(56) 



(1 + f) 3 
R b {r) = 2j-± ^r-<"s 



Af_ _ K 2 \ 
3 24tt i 



4tt 1 



l~2 ^Vy _ 24^ 



(57) 



which are termed as cosmological and black hole horizons. 



II. When F(r) 
apparent horizon 



2(l + /)2 

(\ 2 K 2 \: 

V 3 24 "/ 



(A > jf ) there is only one positive root, which corresponds to a single 



Rbc(r) 



' (1 + /) 3 

(Ai _ K 2 \ 

y 3 24^^ 



(58) 



III. For F(r) > 
apparent horizons. 



2(1+/)!? 

~ K 2 N* 



(A > ~y ) there are no positive roots and consequently there are no 



IV. If A < 



K 

2v / 2~7r 



then also one positive root exists to have a unique apparent horizon. A detailed study of 



apparent horizons in quasi-spherical gravitational collapse with a non-zero A-term can be found in [20] . 



Moreover, solving the dynamical equation (53) we obtain 



R~ 



2V6^F 
V8ttA 2 - K 2 ' 



sinh 



x/8ir\ 2 _K 2 

2_7T 

3 



(t-t s (r)) 



(59) 



where t s (r) is the time of formation of singularity of a particular shell at co-ordinate distance V (i.e., R = 
at t s (r) = 0) . Hence the time difference between the formation of apparent horizon and the singularity formation 
is given by 



t a h{r) -t s (r) 



V8ttA 2 - K 2 



sinh 



V^P" 



2V6^F 



*V 



air 



(One may note that in solving the differential equation (53) it is assumed that A > 



K 



(60) 



However if 



A< 



K 



then 'sinh' is to be replaced by 'sin' in the solution (59)). Thus it is to be observed that heat flow 



and electromagnetic field modify the formation of horizon and also the time difference between the formation 
of apparent horizon and singularity. 

Further, from equation (47) we note that there is dissipation of total energy entrapped inside the surface 
£ due to the heat flow while the electromagnetic field favours entrapped energy within the surface £ . Finally, 
the contribution to the isotropic pressure due to the heat flow and electromagnetic field is in reverse order (or 
in the same order) provided D is positive (or negative). Hence for increasing the electromagnetic field tries 
to diminish the pressure on the boundary which may be balanced by the heat flux term, while for decreasing 
D both of them try to reduce the pressure. So we conclude that electromagnetic field favours formation of 
naked singularity while heat flux may or may not be favourable for formation of black hole. 
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